Potential pitfalls of the PTV concept in dose-to-medium planning optimization

In typical treatment planning of 3D IMRT, the incident energy fluence is optimized to achieve a homogeneous dose distribution to the PTV. The PTV includes the tumour but also healthy tissues that may have a different dose response for the same incident energy fluence, like bony structures included in the PTV (mandibles in head and neck tumours or femoral bones in sarcomas). Dose to medium optimization compensates for this heterogeneous response, leading to a non-homogeneous energy fluence in the PTV and a non-homogeneous dose in the CTV in the presence of geometric errors. We illustrate qualitatively this statement in a cylindrical geometry where the PTV includes a CTV (7 cm diameter) made of water surrounded by ICRU compact bone (1.2 cm thickness); such configuration was chosen to exaggerate the aforementioned effect. Optimization was performed assuming dose equals photon energy fluence times mass energy absorption coefficient. Bone has a 4% lower dose response in a 6 MV flattening filter free spectrum. After optimization either in medium or assuming everything as water composition, the geometry was shifted by 1.2 cm and dose recomputed. As expected, compensating for the under-response of the bone material during optimization in medium leads to an overdosage of the CTV when patient geometric errors are taken into account. Optimization in dose assuming everything as water composition leads to a uniform coverage. Robust optimization or forcing a uniform atomic composition in the PTV margin may resolve this incompatibility between the PTV concept and dose to medium optimization.     

On the choice of the optimal periodic operation for a continuous fermentation process

In this contribution we investigate the impact of the forcing waveform on the productivity of a continuous bioreactor governed by an unstructured, nonlinear kinetic model. The (periodic) forcing is applied on the substrate concentration in the feed. To this end, some alternative waveforms commonly encountered in practice are evaluated and their performance is compared. An analytical/numerical approach is used. The preliminary analytical step is based on the π‐criterion that gives useful information for small amplitudes. The extension to larger amplitudes, when significant improvements are expected, is then performed through a continuation‐optimization procedure. It is found that the choice of the specific waveform has an impact on the performance of the process and there is no unique best forcing for any process condition, but its choice depends on the operating parameters and the forcing amplitude and frequency values. Further, the influence of the waveform functions on the wash‐out conditions are extensively examined. The analysis shows that all the waveforms examined in this work may lead to significant enlargement of the nontrivial regime with respect to a steady state operation. In particular, square‐wave forcing leads in practice to the extinction of the wash‐out conditions for any feed substrate concentration and for a well defined choice of the forcing parameters. © 2010 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

Construction of Selection Indexes with Restrictions

In this paper the problem of unistage selection with inequality constraints is formulated. If the predictor and criterion variables are all normally distributed, this problem can be written as a convex programming problem, with a linear objective function and with linear constraints and a quadratic constraint. Using the duality theory, for convex nonlinear programming it is proved, that its dual problem can be transformed into a convex minimization problem with non-negativity conditions. Good computational methods are known for solving this problem. By the help of the dual problem sufficient conditions for a solution of the original primal problem are derived and illustrated by an example of practical interest.     

Hydraulic acclimation in a Mediterranean oak subjected to permanent throughfall exclusion results in increased stem hydraulic capacitance

Stem water storage capacity and hydraulic capacitance (C_S) play a crucial role in tree survival under drought-stress. To investigate whether C_S adjusts to increasing water deficit, variation in stem water content (StWC) was monitored in vivo for 2 years and related to periodical measurements of tree water potential in Mediterranean Quercus ilex trees subjected either to permanent throughfall exclusion (TE) or to control conditions. Seasonal reductions in StWC were larger in TE trees relative to control ones, resulting in greater seasonal C_S (154 and 80 kg m⁻³ MPa⁻¹, respectively), but only during the first phase of the desorption curve, when predawn water potential was above −1.1 MPa. Below this point, C_S decreased substantially and did not differ between treatments (<20 kg m⁻³ MPa⁻¹). The allometric relationship between tree diameter and sapwood area, measured via electrical resistivity tomography, was not affected by TE. Our results suggest that (a) C_S response to water deficit in the drought-tolerant Q. ilex might be more important to optimize carbon gain during well-hydrated periods than to prevent drought-induced embolism formation during severe drought stress, and (b) enhanced C_S during early summer does not result from proportional increases in sapwood volume, but mostly from increased elastic water.     

On the parametrization of lateral dose profiles in proton radiation therapy

PurposeThe accurate evaluation of the lateral dose profile is an important issue in the field of proton radiation therapy. The beam spread, due to Multiple Coulomb Scattering (MCS), is described by the Molière's theory. To take into account also the contribution of nuclear interactions, modern Treatment Planning Systems (TPSs) generally approximate the dose profiles by a sum of Gaussian functions. In this paper we have compared different parametrizations for the lateral dose profile of protons in water for therapeutical energies: the goal is to improve the performances of the actual treatment planning.MethodsWe have simulated typical dose profiles at the CNAO (Centro Nazionale di Adroterapia Oncologica) beamline with the FLUKA code and validated them with data taken at CNAO considering different energies and depths. We then performed best fits of the lateral dose profiles for different functions using ROOT and MINUIT.ResultsThe accuracy of the best fits was analyzed by evaluating the reduced χ², the number of free parameters of the functions and the calculation time. The best results were obtained with the triple Gaussian and double Gaussian Lorentz–Cauchy functions which have 6 parameters, but good results were also obtained with the so called Gauss–Rutherford function which has only 4 parameters.ConclusionsThe comparison of the studied functions with accurate and validated Monte Carlo calculations and with experimental data from CNAO lead us to propose an original parametrization, the Gauss–Rutherford function, to describe the lateral dose profiles of proton beams.     

基于多尺度协同的长沙市生态网络构建与层级优化

城市化背景下人类与自然环境的矛盾呈现出多尺度、层级化特征，而传统生态网络的构建方式较少考虑不同尺度下生态要素的关系，无法从区域落实到中心城区，难以形成系统性的解决方案。研究在综合梳理各尺度生态网络构建方法的基础上，以长沙市为例，基于形态学空间格局分析（Morphological Spatial Pattern Analysis，MSPA）、景观连通性原理和生态斑块重要性评价识别生态源地，并通过多层级生态阻力面的确定，综合运用最小费用路径（Least-cost path method，LCP）、电路理论、层级传导理论、尺度嵌套等方法对市域、都市区、中心城区的生态网络进行了协同构建和层级优化，最后基于不同尺度生态网络的特点应用并落实到多层级的国土空间规划体系中。研究结果表明：（1）长沙市域生态网络和都市区生态网络具有较好的层级嵌套特征；共识别两尺度生态叠合源地14个、生态叠合廊道15条，主要通过中心城区内的湘江、浏阳河和捞刀河部分河段与外围生态绿圈相衔接，形成"外环内楔"的空间格局。（2）确定市域重要廊道、市域潜在廊道、生态叠合廊道、都市区重要廊道、都市区潜在廊道的核心保护面积共501.14 km²，并提取位于生态廊道核心保护区范围内的生态夹点和生态障碍点，以进一步落实生态保护修复策略。（3）得到具有重要生态连通功能的中心城区生态绿道长度441.2 km，生态修复单元56个，并结合生态阻力值划分为5级进行针对性修复。（4）基于不同尺度生态网络的衔接、嵌套，最终构建"市域总体生态安全格局-都市区城市生态空间发展格局-以城市绿道为基础的中心城区生态修复单元"，并与不同层级的国土空间规划体系相对应。研究结果将为以大城市为中心的跨尺度生态系统修复和生态安全格局构建提供科学参考。     

Effect of shoulder model complexity in upper-body kinematics analysis of the golf swing

The golf swing is a complex full body movement during which the spine and shoulders are highly involved. In order to determine shoulder kinematics during this movement, multibody kinematics optimization (MKO) can be recommended to limit the effect of the soft tissue artifact and to avoid joint dislocations or bone penetration in reconstructed kinematics. Classically, in golf biomechanics research, the shoulder is represented by a 3 degrees-of-freedom model representing the glenohumeral joint. More complex and physiological models are already provided in the scientific literature. Particularly, the model used in this study was a full body model and also described motions of clavicles and scapulae. This study aimed at quantifying the effect of utilizing a more complex and physiological shoulder model when studying the golf swing. Results obtained on 20 golfers showed that a more complex and physiologically-accurate model can more efficiently track experimental markers, which resulted in differences in joint kinematics. Hence, the model with 3 degrees-of-freedom between the humerus and the thorax may be inadequate when combined with MKO and a more physiological model would be beneficial. Finally, results would also be improved through a subject-specific approach for the determination of the segment lengths.     

Development and optimization of a new synthetic process for lorcaserin

A two-step process to synthesize racemic lorcaserin was developed from 2-(4-chlorophenyl)ethanol via formation of bromide or tosylate derivatives. These derivatives were reacted with allylamine in neat conditions to provide pure N-(4-chlorophenethyl)allylammonium chloride. This compound was cyclized in neat conditions using aluminum or zinc chloride to give racemic lorcaserin. After resolution of enantiomers, the wrong enantiomer was racemized and recycled to give new R-lorcaserin.